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ON THE L UNAR AND PLANETAR Y THEORIES. 



BY JOHN N. 8TOCKWELL. 

In my Theory of the Moon's Motion, which was published about a year 
ago, I have called the attention of astronomers to a class of inequalities 
which appear to have been incorrectly calculated by all my predecessors in 
this interesting field of inquiry. The inequalities to which I refer, are of 
great importance in the lunar theory, not only on account of their own mag- 
nitude, but also on account of their modifying influence on some of the oth- 
er inequalities of the moon's motion. They are wholly independent of the 
coordinates of the disturbing body ; and the coefficients of the time in the 
arguments of the equations differ from unity but by quantities of the order 
of the disturbing forces. They arise from the integration of equations of the 
form, 

°^2L. + N 2 dr + m 2 cos(^ — e) = 0; (1) 

in which .Wand i differ from unity but by quantities of the order m 2 , which 
here represents the disturbing function. 
The general integral of equation (1) is 

5r = ^Z^ C08 ^- £ )- ( 2 ) 

For the particular inequality to which I wish to call attention in this pa- 
per, equation (1) becomes 

__T + N*d± — i£m?efcos{nt+a)—28) = ; (3) 

the development being in accordance with Ponte coulant's work on the lu- 
nar theory, and the notation being the same as in my own work ; namely, 
r denotes the radius vector ; e and y, the eccentricity and inclination of the 
moon's orbit; nt, a> and Q } the mean longitude of the moon, the perigee 
and the node, respectively. 

By reason of the perturbations co and 12 are variable, and we have 

dm = —dii = fm'ni. (4) 

The coefficient of nt in the argument of equation (3) therefore becomes 
1-j-fm 2 = i; while iVis equal to uuity minus the motion of the perigee, or 
N=l— $m>. Therefore 

i* — N 2 =-- 6m 2 , (5) 

and equation (3) gives by integration, 
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dl — — f e^cos (nt + to — 23). (6) 

The inequality depending on this argument is, therefore, apparently 
independent of the disturbing force, since the disturbing function becomes 
a factor of both numerator and denominator of the integral expressed by 
equation (2). La Place designates such inequalities as finite. They are, 
however, really indeterminate by equation (2), as we shall now proceed to 
show. It is evident that equation (1) becomes more accurate the smaller 
the disturbing function m 2 becomes, since it is the first term of a converging 
series depending on the ascending powers of m 2 . Equation (1) would there- 
fore be true if m 2 were an infinitesimal ; but equation (2) gives the same 
value of dr whether m 2 be infinitesimal or finite; a result which is in direct 
contravention of the general principle of mechanics that the effect is propor- 
portional to the cause. 

But at this point we are met by the objection of Mr. G. W. Hill, who 
remarks, in the March No. of the Analyst, that the value "m=0 implies 
that we have either an infinitely short month or an infinitely long year" ; 
and is therefore not an admissible hypothesis, since the value of m is really 
finite. 

Now Mr. Hill's remark is correct only in the case in which the sun is the 
disturbing body ; and does not apply to the case of a second satellite or a 
planet. In this latter case, if we put the constant term of the expression of 
the inverse cube of the planet's distance from the earth equal to b-—a" s , a" 
being the mean distance of the planet from the sun ; and also put 
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m*=m"b%- it (7) 

€h 

m" being the planet's mass ; we shall have 

8a) = —33 = f m s «£ and i 2 — iV 2 = 6m 2 ; 
but epuation (6) would remain unchanged. 

In the case of Neptune being the disturbing body, we shall find m 3 = 
soooooV innnns'? an ^ ^is value of the disturbing force would cause the per- 
igee and node of the moon's orbit to move at the rate of about one second 
of arc in a period of 7000 years ; but the substitution of these quantities in 
equation (3) gives the same perturbation as the disturbing force of the sun, 
which is about 500,000,000 times greater. It is therefore evident that eq'n 
(2) fails to give correct results in the case where i differs from unity but by 
quantities of the order m 2 . 

In the lunar theory, however, the disturbing force of the sun is compar- 
atively great in respect to the central force; and as some doubt may exist 
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as to its competency to completely change the character of the elliptical in- 
equalities due to the central force, we shall now show that the application 
of the same method to the planetary theories, in which case the disturbing 
forces are almost infinitesimally small in comparison with the central force, 
leads to results as extravagant and absurd as the perturbations produced by 
analysis (though ascribed to gravita'n), in the theories of the moon's motion. 

For this purpose we shall consider the motion of Mercury as disturbed 
by Neptune. If we designate the coordinates of Mercury and Neptune by 
the same notation which we have employed for the moon and sun, respect- 
ively, in the lunar theory, we shall have in 2 = m'(a 3 -r-a n ) = 0.00000000- 
00114 = the disturbing function. And if we neglect the eccentricity and 
inclination of Neptune's orbit, the perihelion and node of Mercury's orbit 
would have the following motions on the ecliptic, t being reckoned in Juli- 
an years; namely, 

da, == _ dii = + f m 2 ra5 = + 0".0004603 *. 

Now according to Mr. Hill's logic, so long as the perihelion and node 
move at all, equation (2) must give the correct value of the perturbation of 
the radius vector depending on the argument nt-\-co — 242; consequently 

d~ = —&ef cos (nt+a>— 23); 
r 

and this gives for the corresponding perturbation in longitude 

dv = —\ef (= 805") sin (nt+ea— 242). (7) 

Now the pure elliptic motion of Mercury gives rise to the following ine- 
quality in the longitude; namely, 

dv = + ±ef (= 322") sin (nt+co—23); 
and the sum of these two values of dv would be the true value of the equa- 
tion, in order to allow for the perturbation by Neptune. 

Now I find that the greatest perturbation of Mercury by Neptune is that 
depending on the argum't of the evection in the lunar theory, and amounts 
to only 0".0125, while the coefficient of equation (7) amounts to only — 0"- 
.000000193 ; which is altogether more probable than the preceding value. 

Suppose now that we have a second disturbing planet, and call the dis- 
turbing function m' 3 , the motion of the perihelion and node arising in con- 
sequence would be dco = — d3 = fm' 2 n<; and the integral of (3) would be 

8- = — %ef j~ m2 2 +m ^~|cos (nt+io—23) = — fefcos (nt+co—20), 

the same as for a single disturbing planet. It is evident that any number 
of additional disturbing planets might be considered in the same way, and 
give the total perturbation by all the planets the same as that arising from a 
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single one. The perturbations depending on the proposed argument are 
therefore entirely indeterminate by equation (2). 

Now suppose we were to force such an inequality as is given by equation 
(7) into the theory and tables of Mercury ; it is evident that it would give 
rise to several inequalities of sensible magnitude, dependiag on different ar- 
guments; and also that the differences between theory and observation wo'd 
be very nearly equal to the sum of the inequalities arising from the new 
equations thus introduced ; since the existing tables very closely represent 
the motion of Mercury. Assuming the theory including these new equat's 
to be correct, we should endeavor to make the differences between theory 
and observation disappear by applying corrections to the elements of the or- 
bit. In this way we should be able to make the residuals disappear more 
or less completely for a time, since the errors of the elements would partial- 
ly compensate for the errors of the theory and tables. But no amount of 
tinkering of the elements and theory possessing such an inherent source of 
error would suffice to produce tables which would permanently represent 
the motion of Mercury with the precision required by observation. 

Now it appears to me that the lunar theory is in just the condition that 
the theory and tables of Mercury would be in this supposed case. Several 
inequalities of considerable magnitude whieh have no existence in nature 
have been forced in to the theory and tables of the moon ; then corrections of 
the elements have been determined by means of numerous equations of con- 
dition, by which means we have neither correct elements nor correct theo- 
ries ; since equations of condition are powerless to correct for constant or 
systematic sources of error. 

It is proper to remark in this connection, that La Place has stated in book 
II., chapt. V, of the Me'canique Celeste, where the integral of equation (1) 
first appears, that the integral takes a different form from equation (2), in 
the case of N being equal to i. For this supposed case he gives 

3r = -^- 3 cos (Nt-e)-^ sin (Nt-e) (8) 

which is easily proved by differentiation to be correct. Now the last term 
of this equation vanishes at the epoch, when t = ; and we need only con- 
sider the second term in this connection. Since N= 1 — fm 2 , and i = 1 + 
|t» 2 , we shall have N = 0.9958 and i = 1.0126, in the case of the moon 
disturbed by the sun. These numbers are not exactly equal, although very 
nearly so. In the case of Mercury disturbed by Neptune, we have 
JV= 0.99999999999145, and i = 1.00000000002565. 
These numbers approach very nearly to the ratio of equality ; but it may 
very easily be shown that the elements of the orbits should be treated as con- 
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stant in the differential equations, and we should then have, rigorously, 
N = i =■ 1; and the integral of equation (1) would become 

dr = — Jm 2 cos(i< — e). 
Equation (3) would also become 

-J^+dl — fmVcos (nt+at-20) = 0; (9) 

the integral of which would be 

dl = %m?efcos(nt+a>—28). (10) 

This would give for the perturbation of the moon's longitude 

dv = [fmV = 1".2] sin (nt+m—28). (11) 

The perturbation depending on this argument, therefore, amounts to only 
1".2 instead of 111" as determined by Ponte'coulaut, Plana and others. 

From this examination it is apparent that the integral given by equation 
(2) is not applicable to those equations in which the coefficient of t in the 
argument is simply the mean motion of the moon ; and had La Place and 
his successors remembered this circumstance in their calculations, the theory 
of the moon's motion would have been relieved of its most embarrassing 
features; its development would have been confined within narrower limits, 
and its permanent improvement been thereby greatly facilitated. 

Cleveland, Aug. 25, 1882. 

Postscript. — The preceding article was prepared for the November num- 
ber of the Analyst; but as it was not transmitted to the editor until the 
matter for that number had been selected, it was thought best to add the fol- 
lowing as a postscript to that article, in order that the whole difficulty with 
which the lunar theory is embarrassed may be clearly and unmistakably 
traced to its source, 

The equation 

^ +( l_|^_in mV cos (v+w-2S } ~) ^ 

+ lA m ^e'cos(v-a>') J '" ' l } 

given by Plana on page 72, tome II of his The one du Mouvement de la 
Lune, is a particular case of the general equation 

^j[+N 2 y+kcos(it-e) = 0, (13) 

the integral of which is, in the case where N=i, 

y = —jr i cos (its) — -jp sin (it— s). (14) 
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The integral of equ'n (12), neglecting the last term, should therefore be 
du = f^mVcos {v+co— 2J2)-pJ-mV v si » (v+co—2S), (15) 

instead of 

dn = — %efcos (v + co— 2Q\ (16) 

as given by Plana. 

The last term of equation (15) arises from the secular variation of the el- 
meuts of the moon's orbit; while the first term of the second member arises 
from their periodic variations. In order to prove this we shall take the 
values of the secular variation of the elements given by Plana, in tome I, pp. 
96, 97; namely, 

dy = |e 2 r cos (2cp—2Q) 1 . ... „ de = — lefcos (2co—2S) \ n& . 
ydii = |«V sin (2co—2ii) / ' *■ ' edco = ^efsin (2co—2ii) j' { ' 
It is evident that the variations of the elements ought to simultaneously 
vanish at the epoch of the tables ; and as equations (17) and (18) do not sat- 
isfy that condition, we must add a constant to each of these equations in or- 
der to make the variations vanish at a given epoch, when co = co , 8 = £ . 
Equations (17) and (18) will therefore become 

dy = %e\ [cos (2m— 2Q)— cos (2co — 2£ )] 1 , ,-.„, 

ydS = |e 2 r [sin(2«— 2£)-sin (2a> — 2fl )]j' ^ ; 

de — —\ef [cos (2(o— 2S)— cos (2co — 2G j] \ (2Q . 

edco= iey 2 [sm(2o)—2i2)—sm(2w —2ii )]j' <■ ' 

Now the principal term in the value of u is 

u = e cos (v — co ); (21) 

and its variation corresponding to small variations of e and co will be giv- 
en by the equation 

^ w == (-7-)^ e + ( j-)^ = cos ( w — co )de + sin(v — co )edco. (22) 

If we substitute the values of de and edw given by eq'ns (20) we shall find 
du = — ^ey 2 [cos (v — co -|-2a> — 2Q) — cos (v-\- co — 2i? )] 1 ,no\ 

= %ey 2 sin(co — Q — w -fi? )sin (v-\-co — Q — Q ) /' 
In the case of constant elements we have co = co , Q = i2 , and du van- 
ishes. But since co and Q are variable, we have, according to Plana, 

co — co Q +%m\ Si = S — fm 2 «, (24) 

so that co—S2 — ft> -(- j? = ^m?v; and as we are retaining only the terms of 
the order m 2 , we may put 

sin (co — Si — co -\-Si ) — co — Si — co -\-Si = |m 2 v, (25) 

and equation (23) will become 

du = ^-m'ey 2 v sin (v + co — Si — Si — v~\-co — 2Q at the epoch). (26) 
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This is the same as the last term of equation (15), and proves the first p't 
of the proposition. The second part may be proved in a similar manner 
by taking the variation of the second term of the value of u and substitu- 
ting the periodic variation of the elements. 

We thus see that the variation of the elements leads to the same results 
as are derived from the variation of the coordinates, when the proper con- 
stants are added to the integrals of these variations. 

The integral of equation (12) depending on cos (v—co') is 

3u = — -J-f m 2 (a-^-a') e'cos (v — co') — ^-fm 2 (a-ha')e' v sin (v — co'), (27) 
instead of Su — -f (a-i-a') e'cos (v — co'), as given by Plana. 

The last term of equation (27) arises from the secular variation of e and 
co, depending on the distance between the perigees of the sun and moon. To 
prove this we shall observe that the differential equations of e and co contain 
the terms 

= — Ifm 2 -^ e'sin (to— at') \ 

f v a y. (28) 

e— = — ffwr— - e'cos (co — co') 
dv a' J 

If we integrate these equations in the same way that Plana has done, and 
also add a constant to the integrals so that de and dco may simultaneously 
vanish at the epoch, we shall find, 

de = f (a-i-a') e'[cos (co— a>')— cos (to — a/ )] 1 _ ^ 

edco = — f (a-^a') e'[sin (co — co') — sin (w — w' )] J ' 

If these values be substituted in equation (22), we shall find, after making 
the necessary reductions, 

du = — \^m 2 (a-r-a') e' v sin (v — co'), (30) 

as in the last term of equation (27). 

The other term of equation (27) may be found in like manner by substi- 
tuting the proper periodic variations of e and co in equation (22). 

The two cases which we have examined are the most important ones 
among the equations of the moon's longitude; but there are two terms of 
considerable importance in the equations of the moon's latitude, which we 
shall now examine. 

The first of these terms of the latitude is given by Plana in tome II, page 
118, of his work, as follows: 

^ir+(l +|m 2 )&- > 2 eV sin (v— 2co+ Q) = 0. (31) 

The integral of this equation is 

ds = ^j-f-mVj' sin (v — 2co-\-Q) — l-fmWyv cos (v — 2co-\-S); (32) 
instead of ds = 4-fe 2 j-sin (» — 2co-\-Q), as given by Plana. 
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The last term of equation (32) is produced by the secular variation of the 
node and inclination of the moon's orbit, as may be shown in the following 
manner:— Since the principal term of the latitude is given by the equation 

s = rsin(v—£ ), (33) 

we shall have 

8s = Bm(v—S )8 r —oos(v-S )rdS; (34) 

and if we substitute the values of dy and yd®, which are given by eq. (19) 
we shall find 

ds = feYsin (<o — Q Q — io-\-ii) cos (y — a> — co-j-ti) 1 ,„,-x 
= — J^mVj- v cos (v— 2o) + Q) }' 

the same as the last term of equation (32). 

The second term of the latitude which seems to be entirely erroneous is 
the one arising from the oblateness of the earth. La Place first computed 
this inequality ; and I am especially desirous of calling attention to it, since 
the error is common to every theory of the moon's motion which I have ex- 
amined, including my own. In this examination I shall first follow the 
method given by La Place in the Me'oanique Celeste, making use of the 
translation by Bowditch, since it is more convenient for reference than the 
original work. I shall then show how the same results may be obtained by 
the method given in my own work. 

If, in the function [5347], Me'o. Ce'l., we substitute for u-i-h 2 its value 
1-i-a 2 , and also put for abridgement 

m = {ap—laf)—, (36) 

it will become 

2m sin £ cos £ sin fv-\-(g 2 — l)Hsmfv, (37) 

£ denoting the obliquity of the ecliptic to the equator. 
But in [5350] La Place gives 

H = s sin e cos £ : (38) 

9* — 1 

consequently the function (37) is identically equal to nothing. 

Now La Place remarks that the oblateness of the earth adds the function 

(37) to the differential equation of the latitude, by which means it becomes 

dNh 

dv 

the integral of which he gives as 

ds = sin £ cos £ sin fv. (40) 

if 

In other words, if nothing be added to the differential equation of the lat- 
itude, its integral will be increased by the value of ds in equation (40) ; — 
a result which is exceedingly interesting as well as mysterious to a person 
unacquainted with the wonderful powers of the calculus. 



y +&s+[2m sin £ cos e + {g 2 —l)HJ sinfv = 0; (39) 
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But the correct value of the integral of equation (39) is 

8s = — J[2ra sine cos e+C^ 2 — l)H~]amfo 1 ,.-.^ 

-f|[2msin ecoss-)-^ 2 — l)-ff] v cosfv J ' ' 

since /= 1. 

The function H is therefore not correctly given by equation (38), and the 
function (37) is not equal to zero. 

In equation (41) the term — |msin s cos esin/» gives the direct effect of 
the earth's oblateness on the moon's latitude; while the term (g 2 -l)H sinfo 
is its indirect action transmitted to the moon by means of the sun but de- 
creased in the ratio of g z — 1 to unity. Now Jmsin scos e = 0".01654 
= H, and (# 2 — l)H= 0".000133. Whence it follows that the effect of the 
earth's oblateness on the moon's latitude is insensible. 

The same result may be obtained by the method empl'd in my own work. 
For the same force which gives rise to the first term of the function (37), 

produces in the value of " , eq'n (263) of my Theory of the Moon's Mo- 

OX/ 

tion, the following terms 

— ^— = mn sin e cos e [J cos nt — nt sin nt~] ; (42) 

and this gives by integration 

o d = m sin s cos e [ — Jsin nt-\-nt cos ni], (43) 

which is the same as the corresponding terms of ds in equation (41), since 
v = nt, when we neglect the eccentricity of the moon's orbit. 

It is easy to prove that equation (43) is correct, by means of the variation 
of the elements. For if we add a constant to the secular terms of these va- 
riations, which are given by equations (630) and (631) of my work, so that 
they may simultaneously vanish at the epoch, they will become 

j8Q = m sin e cos s [Jsin (2nt—@) — (n-z-a') (sin Q— sin X? ) \ .... 

by = m sin s cos e [^cos (2nt— ■S)-)-(n-s-a / ) (cos Q— cos S Q ) j ' ' > 
And if we substitute these values in the equation 

8d = sin (nt— 8 )o?— cos (nt—S )ydS, (45) 

the secular terms will produce the following term, 

dd = 2 — r sin e cose sin \(Q — i? )cosn£. (46) 

But since S2-Q Q —a% we have sin \(Q-Q Q ) = \a% and eq. (46) becomes 
dd = m sin s cos e nt cos nt; (47) 

which is the same as the last term of (43). 

If we now substitute the periodic terms of equations (44) in equation (45) 
we get the first term of equation (43). 

Finally La Place remarks, Me'o. Ce'L, [5398], that the inequality of the 
moon's motion in latitude arising from the oblateness of the earth is only 
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the reaction of the nutation of the earth's axis, discovered by Bradley ; — a 
statement which seems strangely at variance with mechanical principles, 
since the one has a period of only twenty-seven days, while the other has a 
period of nearly nineteen years. 

It is evident that this change in the equation of the moon's latitude affects 
the corresponding equation of the longitude ; and I find that it changes its 

val ue from dv = 6 — =- r sin e cos e sin ii, to dv = -V£ — r r sin s cos e sin ii. 

a' ' ' 2 a' ' 

The equations of the moon's motion arising from the oblateness of the earth 

are, therefore 

dv = 4".814 sin Q, dO = — 0".0165 sin nt, (48) 

instead of 

dv = 4".444 sin 12, dd = — 8".226 sin nt, 

previously found. 

But there is one point in Mr. Hill's "Review" which seems to be well 
made out, and for which I wish to thus publicly tender him my grateful 
acknowledgements. And that is in relation to the reduction of the moon's 
longitude from the plane of its orbit to the ecliptic, and which amounts to 
the correction which I have now applied to the equation of the longitude 
arising from the oblateness of the earth. 

I have thus shown with considerable detail, not only that the lunar the- 
ories of my predecessors are erroneous by terms of the third order, but have 
also shown just how the errors were introduced into the several theories. 
The errors were first committed by La Place; and his methods of integra- 
tion were blindly followed by his successors until his results came to be 
accepted by astronomers as the legitimate effects of the law of universal 
gravitation. Had La Place been a less skilful mathematician it is not prob- 
able that he would have made such mistakes; since he would have been 
content to consider one power of the disturbing force at a time, and he wo'd 
then have been obliged to regard the elements as constant in the integrates ; 
and then, in all the cases which I have pointed out, he would have found 
zero for a divisor and infinity for the perturbations. And such results 
could not have failed to indicate their origin in the improper application of 
a general formula of integration to a special case in which such formula fails 
to give correct results. I therefore confidently believe that existing theories 
when corrected for the sources of error which I have pointed out; together 
with a careful revision of the elem'tsof the lunar orbit, will give the moon's 
place in the heavens with all the precision required by observation, without 
carrying the approximat'n to terms of a higher order than have been already 
calculated. 

Cleveland, Oct. 10, 1882. 



